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ABSTRACT 

ISO-SWS spectroscopy of the WR galaxy NGC 5253 is presented, and analysed to pro- 
>• vide estimates of its hot young star population. Our approach differs from previous 

investigations in that we are able to distinguish between the regions in which different 
infrared fine-structure lines form, using complementary ground-based observations. 
The high excitation nebular [S iv] emission is formed in a very compact region, which 
we attribute to the central super-star-nucleus, and lower excitation [Neil] nebular 
emission originates in the galactic core. We use photo-ionization modelling coupled 
with the latest theoretical O-star flux distributions to derive effective stellar tempera- 
tures and ionization parameters of T e ff > 38kK, log Q ~ 8.25 for the compact nucleus, 
with T e ff ^35kK, log Q <8 for the larger core. Results are supported by more sophis- 
ticated calculations using evolutionary synthesis models. We assess the contribution 
q ■ that Wolf-Rayet stars may make to highly ionized nebular lines (e.g. [Oiv]). 

?— i ' From our Bra flux, the 2" nucleus contains the equivalent of approximately 1 000 

^ \ 07 V star equivalents and the starburst there is 2-3 Myr old; the 20" core contains 

■ about 2 500 07 V star equivalents, with a representative age of ~5Myr. The Lyman 

ionizing flux of the nucleus is equivalent to the 30 Doradus region. These quantities 
are in good agreement with the observed mid-IR dust luminosity of 7.8 x 10 8 Lq. Since 

■ this structure of hot clusters embedded in cooler emission may be common in dwarf 

starbursts, observing a galaxy solely with a large aperture may result in confusion. 
Neglecting the spatial distribution of nebular emission in NGC 5253, implies 'global' 
stellar temperatures (or ages) of 36kK (4.8 Myr) and 39kK (2.9 or 4.4 Myr) from the 
observed [Nem/ii] and [Siv/iii] line ratios, assuming log Q=8. 

Key words: galaxies: individual (NGC 5253) — galaxies: starburst - galaxies: stellar 
content - infrared: galaxies - stars: Wolf-Rayet 
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1 INTRODUCTION 

Wolf-Rayet (WR) galaxies are a subset of blue emission-line 
galaxies whose spectra show the signature of large numbers 
of WR stars (Vacca & Conti 1992). These galaxies are among 
the youngest starburst galaxies and contain a large popula- 
tion of very massive stars with ages in the range 1-10 xlO 6 
yr. Mid infra-red (IR) observations are well suited to the 
study of hot, massive stars since this spectral region contains 
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many fine-structure nebular lines which depend very sensi- 
tively on stellar content. Until recently, use of such diagnos- 
tics was severely hindered because of the low transparency 
of our atmosphere at mid-IR wavelengths. The launch of the 
Infrared Space Observatory (ISO) has opened up this new 
window, which we are exploiting through a Guest Observer 
program (WRGALXS; P.I. A.J.Willis) with the ISO Short 
Wavelength Spectrometer (SWS). 

NGC 5253 is one of the closest WR galaxies, at a dis- 
tance of only 4.1 Mpc (Saha et al. 1995) and represents one of 
the targets of our Guest Observer program. NGC 5253 con- 
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Figure 1. ISO— SWS apertures superimposed on the Calzetti et 
al. (1997) HST/WFPC2 Ha image of NGC 5253, together with 
the 2 cm VLA radio continuum contour map from Beck et al. 
(1996). The ISO aperture is 14x20 arcsec for A<12/mi, and 14x27 
arcsec at longer wavelengths. The spacecraft Y- and Z-axes, indi- 
cated here, are discussed in the text. Declination (—31 38m) and 
right ascension (13h 39m) are in J2000.0 coordinates. 

tains several very young super-star clusters, no older than a 
few million years (Walsh & Roy 1989; Beck et al. 1996). It 
has been the focus of many recent studies (e.g. Gorjian 1996; 
Calzetti et al. 1997; Turner, Beck & Hurt 1997), which have 
shown that the extinction is very high and patchy at ultravi- 
olet and visual wavelengths. ISO observations of NGC 5253 
are particularly important since the interstellar extinction at 
mid-IR wavelengths is low and because earlier spectral ob- 
servations found evidence for a remarkably hot and young 
stellar population (Aitken et al. 1982, Lutz et al. 1996; Beck 
et al. 1996). 

In Section ^ we present our new ISO-SWS spec- 
troscopy of NGC 5253 together with complementary ground- 
based observations, which are interpreted in Section ^[ Our 
technique for photo-ionization modelling, based largely on 
CLOUDY (Ferland 1996) is discussed in Section ^, including 
assumed elemental abundances and theoretical flux distri- 
butions for O stars. In Section ^ we discuss the results of 
our photo-ionization modelling. The influence from WN and 
WC-type Wolf-Rayet stars is considered in Section ^| while 
results obtained using evolutionary synthesis models are pre- 
sented in Section Our present results are compared with 
those from the literature in Section ^ and our conclusions 
are reached in Section fol 



2 OBSERVATIONS AND REDUCTION 

Our new infrared data on NGC 5253 were obtained as part 
of a Guest Observer programme with the Short Wavelength 



Spectrograph (SWS; de Graauw et al. 1996) onboard the 
ESA Infrared Space Observatory (Kessler et al. 1996) in 
revolution 418, 7 January 1997. 

The SWS AOT6 observing mode was used to achieve 
full grating resolution, A/ A A ~ 1300 — 2500 and the contin- 
uous wavelength coverage was 2.99-27.65 /im. Data of detec- 
tor "bands" 1 and 2 respectively cover wavelengths of 2.99- 
4.08 pm using 12 In:Sb detectors and 4.00-12.05 [im with 12 
Si:Ga detectors, employing entrance slits that give an effec- 
tive aperture area of 14 x 20 arcsec on the sky. Band 3 data 
cover 12.0-27.65 /J,m using 12 Si:As detectors, with sky cov- 
erage of 14 x 27 arcsec. The total integration time was set 
to 9140 seconds to allow one complete scan over the wave- 
lengths selected within each "AOT band", defined by the 
permissible combinations of detector band, aperture, and 
spectral order (cf. de Graauw et al. 1996). The observation 
time includes dark current measurements, and a monitor 
of photometric drift for the detectors of bands 2-3. A drift 
measurement is not normally made for the relatively stable 
In:Sb detectors. 

The SWS data were processed through the main stages 
of raw signal conversion (bits to pV/sec), glitch recogni- 
tion, wavelength calibration, and dark current subtraction 
by following the standard chain of pipeline software and cal- 
ibration tables complying with version 6.0 or later. Final 
photometric calibrations performed in the SWS Interactive 
Analysis environment rely on relative spectral response func- 
tions and absolute photometric response factors that were 
conditional for further pipeline development. These are gen- 
erally consistent with calibrations found in pipeline version 
7.0. Prior to co-addition of each detector's output, data were 
statistically filtered using the results of glitch recognition al- 
gorithms, rejecting data points whose corresponding signal 
ramps were found to have more than two departures from 
linearity after AC correction. Detector 27 was judged to be- 
have photometrically unstable in this observation, and all 
output from it was rejected. Remaining detector output was 
rectified to the mean flux density within each resolution el- 
ements, 3-cr points were excluded, and finally re-binned to 
a spectral resolution that is somewhat higher than instru- 
mental. 

General aspects of the wavelength and photometric cali- 
brations as of revolution ~100 (after initial verification of in- 
strument performance) are described, respectively, by Valen- 
tijn et al. (1996) and Schaeidt et al. (1996). More contem- 
porary performance details are in preparation for the ISO 
Post Mission Archive by the SWS instrument team. 

2.1 Aperture orientation 

The input coordinates coinciding with aperture centre were 
q = 13:39:55.7, 5 = -31:38:29.0 (J2000.0). Spacecraft point- 
ing uncertainties are alleged to be around 2 arcsec at the 
time of this pointing, following recalibration of the star- 
tracker CCD some 50 revolutions prior. Pointing reconstruc- 
tion software in pipeline version 7.0, accounting also for later 
changes in solar ephemerides updating procedures, indicates 
offsets, both along dispersed (corresponding to the space- 
craft Z axis) and orthogonal to the dispersed direction of 
SWS (Y axis) to be within the quoted uncertainty. The cen- 
tre of our apertures is ~2" from the ultraviolet source UV3 
of Kobulnicky et al. (1997), but is sufficiently large to in- 
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elude the strong nebular Ha and Bra emitting regions of 
NGC 5253, ~5" from the centre of our aperture (see Kobul- 
nicky et al. 1997; Davies, Sugai & Ward 1998). The two SWS 
aperture projections are overlaid on a Hubble Space Tele- 
scope (HST) WFPC2 Ha image of NGC 5253 from Calzetti 
et al. (1997) in Fig. jj], together with a 2 cm VLA contour 
map from Beck et al. (1996). 

It is important to note that photometric responses are 
not flat over the spatial extent of the area covered by each 
detector block. This is most important in the direction or- 
thogonal to the dispersion axis (along the spacecraft Y axis; 
see Fig. |l|), where an offset of 5 arcsec from the position of 
peak responsivity results in a loss of 20-25% flux of com- 
bined detector output from a point source. The loss within 
2 arcsec of the aperture centre is less than 5% for bands 
2-3, which exhibit similar beam profiles. This is shown in 
rough form in Schaeidt et al. (1996), but detailed profiles to 
aperture edge and beyond are in analysis at this time. For 
band 1, the loss is negligible within ~2.5 arcsec of the aper- 
ture centre, but a sharper gradient in responsivity beyond 
this again results in a loss of 20-25% at ~5 arcsec away. 
If 30% of Bra line emission does originate from the central 
nucleus of the NGC 5253 (cf. Davies et al. 1998), then the 
offset along the minor axis of the aperture may mean that 
the energy measured in the line is underestimated by ~10%. 
For more extended line emission (see below), a loss estimate 
is more difficult to arrive at without a surface brightness dis- 
tribution and knowledge of the photometric responses to the 
aperture edges and beyond. Overall, photometric uncertain- 
ties are estimated to be ~5% in band 1, 10-15% in band 
2, and 15-25% in band 3. No large discrepancies between 
levels of adjoining AOT bands were noticed. 



2.2 Infrared nebular emission line fluxes 

We present the full SWS spectrum of NGC 5253 in Fig. |[ 
Several sharp nebular emission lines are apparent, su- 
perimposed on a rising 'continuum' flux distribution at- 



Table 1. Nebular emission lines observed in our ISO— SWS 
(AOT6) spectrum of NGC 5253, and ground-based measure- 
ments. Our ISO fluxes are after correction for the v aria ble profile 
across the beam in the spacecraft Y-axis (see Sect. j2,l|). 



tributed to dust emission (see Section 3.2). The most pro- 
nounced nebular lines are due to ground-state fine structure 
lines of [Siv]10.51/xm, [Ne ill] 15.55pm and [S lll]18.71/xm; 
weaker nebular emissions are also seen in [Nell]12.81/im and 
[Ar lll]8.99pm. We have used the emission line fitting routine 
elf in DIPSO (Howarth et al. 1995) to measure line fluxes 
using an eye-estimated local continuum flux, which are pre- 
sented in Table Since the ISO aperture is offset from the 
central nucleus along the minor axis, we have included a 
10-15% correction for light loss. From the Table, there are 
clearly additional uncertainties - the measured [S iv] flux 
from the small aperture ground-based observation is greater 
than that from the large aperture ISO data set. 

Individual nebular line profiles are shown in Fig |^, al- 
though these are unresolved at the SWS AOT6 resolving 
power. Emission lines are red-shifted by Ai;~+200 km s , 
in reasonable agreement with optical measurements of the 
galactic recession velocity (+416 kms -1 from de Valcouleurs 
et al. 1993). The observed Av was then used to search for 
weaker emission lines allowing us to confirm the presence 
of Br a 4.051pm and (possibly) Pfa 7.458pm (see Fig §). 
We include the spectral region around [Arll]6.99pm in our 
band 2C SWS dataset, although the positive identification 
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and measured flux of this feature is uncertain, due to its 
location at the edge of the 2B and 2C bands and low Av. 

We have also collected IR flux measurements for 
NGC 5253 from the literature which, are also included in 
Table [j]. The hydrogen Brackett line data are taken from 
Kawara, Nishida & Phillips (1989) and Davies et al. (1998) 
and the metal lines from Beck et al. (1996). In general the 
ISO-SWS and ground-based results agree reasonably well, 
except for the [Neii]12.81/im line (see Section^). Genzel et 
al. (1998) have recently presented additional [S ill] and [Sill] 
mid-IR line flux ratios obtained with band 4 of ISO/SWS, 
obtained at a location ~2.4 arcsec from our data set. We 
have converted these into absolute flux measurements us- 
ing our [Nell] measurement, although this introduces an 
additional uncertainty since the SWS band 4 aperture is 
20x33 arcsec. In theory we should be able to use the ob- 
served ratio of [S ill] (33.5pm/18.7pm) of ~1.6 to measure 
the (emissivity weighted) mean electron density of the [S ill] 
emitting region. However, these lines were obtained through 
different sized apertures, and the observed ratio is close to 
the low density limit, so applying simple nebular theory, we 
find n e <10 2_3 cm~ 3 . For comparison, Walsh & Roy (1989) 
obtained n,; = 10 2 ' 46 cm -3 for the central ~ 8 X 4 arcsec of 
NGC 5253 from optical spectroscopy. 



3 INTERPRETATION OF GROUND AND 
SPACE-BASED IR OBSERVATIONS 

Before we undertake a detailed analysis of our ISO observa- 
tions we first compare our ground and space-based mid-IR 
line fluxes. We also estimate the massive stellar content of 
NGC 5253 from the IR dust continuum. There are potential 
dangers in forming line ratios when the source is a galaxy, 
which can be expected to contain many regions of different 
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Figure 2. Flux calibrated (in Jansky) spectrum of NGC 5253 obtained with ISO/SWS, including a black body fit (dotted line) to the 
observed dust continuum. The broad feature between 22-27/^m appears to be due to solid state silicate emission, representing the first 
such detection in a starburst galaxy. 



characteristics, and especially if it is a galaxy like NGC 5253 
which is known for marked spatial variations in stellar age, 
dust content and other features. 



3.1 A Hot Nucleus in a Cooler Extended Core 

The infrared line fluxes listed in Table |l] include both the 
ISO results and, for the [Siv], [Nell], and [Arm] lines, the 
fluxes observed from the ground with the Irshell spectrom- 
eter (Beck et al. 1996) which had a 1.6" beam; mapping 
showed that the [S iv] and [Ar ill] detected by Irshell came 
from a region about 2" in diameter. Similar dimensions were 
found by Davies et al. (1998) for the dominant Br7 source. 
At the distance of NGC 5253, this corresponds to a (maxi- 
mum) size of 40 parsec. The most striking thing about Ta- 
ble |l| is the close agreement of the [S iv] 10.5/xm fluxes ob- 
served with the large ISO beam and the 1.6" Irshell beam. 
This result argues that the [S iv] emission is concentrated 
in the inner 2-4" of the galaxy. However, we cannot assume 
that all the line fluxes are similarly concentrated; another 
feature of Table |l| is the order of magnitude excess of [Ne n] 
12.8/im emission in the ISO observations compared to the 
ground-based, which is most readily explained if the [Nell] 
is extended over an area comparable to the ISO beam so the 
small Irshell beam did not measure the total. 

In Section ^| we shall show that [S iv] emission increases 
with increasing stellar temperatures and ionization param- 
eters, with the reverse found for [Nell]. Therefore the most 



natural explanation is that the [S iv] region, the central 2" 
mapped by Irshell, contains very hot stars at high ionization 
parameter, while the [Nell] region contains cooler stars at 
lower Q-value. Other lines from high ionization stages will 
be likely to originate in the [S iv] region and those from lower 
ionization in the [Ne n] region. 

This is a far-reaching assumption, but the picture of 
NGC 5253 that has been established from high spatial reso- 
lution observations supports it. In the radio continuum the 
thermal emission due to young stars is dominated by a cen- 
tral source of less than 2" which contains more than half 
the radio flux from the area of the ISO beam (Turner et al. 
1998). This radio source is extremely young, which has led 
it to be identified with the 'central star cluster' [S iv] source 
(Beck et al. 1996) and with the visually obscured cluster 
NGC 5253-5 of Calzetti et al. (1997) which is thought to be 
no more than 2.5 Myr old, resides in the center of the star- 
burst nucleus, and is very compact. The youth of the cluster 
is seen in the lack of supernovae remnants and the presence 
of very massive and shortlived stars which are hot enough 
to excite strong [S iv] emission. Calzetti et al. distinguish 
between the starburst nucleus, which contains the central 
cluster, and the galactic core, which is older, produces stars 
at a rate about 1/10 that of the nucleus, is UV rich, and cov- 
ers an area comparable to the ISO beam. If the [S iv] source 
is associated with the starburst nucleus, the [Ne n] emission 
may be presumed to come from the larger core. The halo of 
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Figure 3. Principal nebular emission lines observed in the ISO— 
SWS spectrum of NGC 5253, red-shifted by Av — (-200 km s -1 . 
The identification of [Aril] and Pf a are uncertain. 



NGC 5253 was not observed by ISO and is not discussed in 
this work. 



3.2 Dust Emission and the Number of O stars in 
NGC 5253 

The SWS spectrum of NGC 5253 in Fig | shows that the 
nebular emission lines are superimposed on an underlying 
energy distribution which has negligible flux below 4 /im, 
peaks at about 10 Jy at 24 fim and then turns over to longer 
wavelengths. This energy distribution clearly reflects emis- 
sion from relatively warm dust. The broad emission between 
22.5-27/im suggests possible solid state silicate emission, a 
sensitive tracer of mass-loss history. This is common in O- 
rich environments of the Galaxy but is not known in extra- 
galactic regions, like 30 Doradus, where the dust may be in- 
sufficiently shielded. Therefore, NGC 5253 appears to be the 
first example of silicate emission in a starburst galaxy. Us- 
ing blackbody energy distributions normalised to the SWS 
data at about 21/im we find a best fit to the observed 
data for a temperature of 200±10 K, as shown in Fig. ^| 
From this, the total integrated dust emission at the earth 
is 1.6xlO -12 Js -1 m -2 . Adopting a distance of 4.1 Mpc to 
NGC 5253 (Saha et al. 1995), this translates into a total 
dust emission luminosity of 2.9xl0 35 Js -1 , or Ltr=7.8x10 8 
L.0. (For comparison, Lutz et al. (1996) report Lir =8.0xl0 8 



Table 2. Stellar parameters and ionizing fluxes of main-sequence 
O star flux distributions from Schaerer & de Koter (1997), for 0.2 
Zq . Approximate spectral types follow the temperature calibra- 
tions of Vacca et al. (1996, WDV) and Crowther (1997, PAC) 



Mass 


Age 


T e ff 


logL* 


logQH 


logQHd lo 


g<?HeII 


Sp. Type 




Mq 


Myr 


kK 


Lq 


phs -1 


ph s — 1 


ph s — 1 


WDV PAC 


A2 


20 


3.6 


33.3 


4.77 


47.84 


46.08 


0.00 


B0V 08.5V 


132 


25 


2.6 


36.3 


5.00 


48.43 


47.42 


44.23 


09V 08V 


C2 


10 


1.5 


41.7 


5.45 


49.11 


48.53 


46.04 


07V 06V 


D2 


60 


0.8 


46.1 


5.77 


49.51 


49.02 


46.45 


05V 04V 


E2 


85 


0.7 


48.5 


6.04 


49.82 


49.36 


46.59 


04V 03V 


F2 


120 


0.6 


50.3 


6.26 


50.07 


49.62 


47.00 


03V 03V 



Lq from ISO- SWS spectroscopy, while Stevens & Strickland 
(1998) derive 6.9x10 s Lq using IRAS flux measurements 
plus our adopted distance). 

Before we estimate the number of hot stars exciting 
this dust, let us recall the concept of an "equivalent" O star 
(Vacca 1994). This hypothetical object was introduced so 
as to ease the transformation between what is measured by 
the Lyman continuum flux and "the number of O stars" and 
provides a useful means of comparing the energetics of star 
formation regions in galaxies and other sources. The number 
of Lyman continuum photons from an O star "equivalent" 
is 10 49 Ly photons -1 , close to that of an 07 V star. Vacca 
(1994) calculates the correction in transforming from the 
number of "equivalent" O stars to the "number of O stars" . 
This depends on the initial mass function (IMF) slope, the 
metal abundance, and the upper mass cut-off. Roughly half 
of the ionizing flux comes from stars hotter and the other 
half from stars cooler than the O star "equivalent', thus 
the correction factor is typically less than a factor two. It 
should be kept in mind that one might be measuring the 
Lyman continuum flux and determining the number of O 
star "equivalents" , and estimating the "number of O stars" 
but there may only be stars of later spectral type than 07 V 
if the starburst is somewhat evolved. In fact, we will see this 
case below in the "core" of NGC 5253. 

The Lyman continuum flux of a region of hot stars is 
usually found from the extinction-free radio continuum, but 
in NGC 5253 this is partly optically thick and there is in- 
sufficient spatial resolution to determine which can be as- 
sociated with the starburst nucleus and which with the 
core. We therefore assume that the Hn region is dust- 
free, and optically thick to Lyman continuum radiation 
(Case B ionization-bounded) . We adopt A r e = 10 2 ' 5 cm -3 and 
T e =ll,700K for the core (Walsh & Roy 1989) for the de- 
termination of hydrogen intensity ratios (Storey & Hummer 
1995). Our extinction-corrected ISO Bra flux implies a total 
ionizing flux of 3.45xl0° 2 Ly photons -1 . This supports the 
previous determination of 4xl0 52 Ly photons -1 from Ha 
by Martin & Kennicutt (1995) though Davies et al. (1998) 
recently used Br7 observations plus an aperture of 30 arcsec 
to derive a factor of two higher Lyman ionizing flux (scaled 
to our assumed distance) . This discrepancy would be larger 
still (factor of 2.7) adopting a consistent IR extinction. 

Davies et al. revealed a dominant nebular source, spa- 
tially coincident with the visual Ha emission peak, contain- 
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ing 30% of the Br7 emission within 2" of the central nucleus. 
(A similar fraction was measured by Martin & Kennicutt 
(1995) from Ha observations.) We therefore conclude that 
the ionizing flux of the nucleus, as derived from our mea- 
sured Bra flux, is 1.0x10 s2 Ly photons -1 - comparable to 
30 Doradus (Kennicutt 1984) - and that the ionizing flux in 
the core is 2.5xl0 52 Ly photons -1 . The total ionization is 
equal to 3 500 equivalent 07 V stars. 

Using the recent O star calibration of Vacca, Garmany 
& Shull (1996), we find a luminosity of ~2.5xl0 5 L for 
an 07 V star. With the typical assumption that the dust 
luminosity is governed by the radiation from the hot stars, 
our derived Lir is equivalent to about 3 100 07 V stars. This 
number is in good agreement with that for the O star "equiv- 
alents" inferred from the Lyman continuum measures. 



4 PHOTO-IONIZATION CALCULATIONS 

The ISO data on NGC 5253 include the important diagnos- 
tic lines of sulphur, argon and neon, from whose ratios it 
should be possible to deduce the extreme-ultraviolet (EUV) 
spectrum of the input ionising radiation, and thus the tem- 
perature of the exciting star (in a nebula excited by a single 
star) or the effective temperature and thus the age and mass 
function of an exciting cluster of stars. In the following sec- 
tion we will use results from a grid of model Hll regions 
produced with the photo-ionization code CLOUDY (Ferland 
1996). First, we describe our technique, choice of elemen- 
tal abundances and O star energy distributions, since these 
quantities will have a major effect on the results. 



4.1 Description of the calculations 

We constructed photo-ionization models using CLOUDY 
(v90.04) as described in Ferland (1996) and Ferland et al. 
(1998). Comparisons with other photo-ionization codes are 
provided by Ferland et al. (1995). The nebula are repre- 
sented by a sphere of uniform gas density, n, and filling 
factor, e, with a small central cavity, which is ionized and 
heated solely by the UV radiation of a single central star. 
Nebular fluxes are predicted, given input abundances, flux 
distributions and physical parameters, most important of 
which is the ionization parameter. This may be written as 
Q = Qu/(iirRgn), where Qh is the number of ionizing pho- 
tons below the H-Lyman edge at 912A, and Rs is the radius 
of the Stromgren sphere. (An alternative definition of the 
ionization parameter is U = Q/c.) 

For a given energy distribution of the ionizing radiation 
field, any combination of parameters which keeps Qnne 2 
constant will result in an identical ionization structure of 
the gas (see Stasifiska & Leitherer 1996). For each flux dis- 
tribution we have followed Stasifiska & Schaerer (1997) by 
computing a series of models with hydrogen density n=10 
cm -3 and e—1, in which the stellar luminosity has been mul- 
tiplied by a factor of 10 5 , 10 4 , 10 -3 . These models thus 
correspond to varying the ionization parameter log Q = 6 to 
10. This covers the range of ionization parameters expected 
in astophysically relevant conditions as found by Stasifiska 
& Leitherer (1996). 



4.2 Abundance Effects 

The elemental abundances in the ionized gas can strongly 
affect the cooling of the nebula and thus the observed line ra- 
tios. NGC 5253 is seen from optical spectroscopy by Walsh 
& Roy (1989) to be metal-poor, with 12+log(0/H)=8.17 
derived for the region centred on the major Ha (and Br7) 
emission. More recently, Kobulnicky et al. (1997) have de- 
rived elemental abundances of N, O, Ne, S and Si in sev- 
eral regions of the core H II region from HST spectroscopy 
and confirm a depleted metal content. They supported ear- 
lier evidence for a region of enhanced nitrogen abundance, 
attributed to enriched material from WN stars. The abun- 
dance of sulphur is 1/4-1/2 times that of the Orion nebula 
(Baldwin et al. 1991). For our analysis we use abundances 
derived for region Hll-1 by Kobulnicky et al. (1997), and 
1/3 that of the Orion nebula (or approximately 1/4 solar 
metallicity) for elements Kobulnicky et al. were unable to 
measure. 

Test calculations were performed including varying of 
dust grains abundances and compositions (ISM, Orion). 
Overall the influence of dust was not substantial, with the 
greatest differences obtained in high luminosity models, such 
that the ionization balance and nebular temperature in- 
crease. We found that the optimum consistency between 
different IR diagnostics was obtained for models including 
ISM grains (silicate and graphite). Consequently, all models 
considered in this work assume this grain mixture. 

4.3 Flux distributions for O-type stars 

Most of the recent studies on analysis of Hll regions use 
the Kurucz (1991) plane-parallel LTE model atmospheres. 
This is because these models include line blanketing and are 
available for a wide range of stellar temperatures and metal- 
licities. However, with regard to O-type stars, Kurucz mod- 
els neglect both non-LTE effects and the influence of stellar 
winds, with important consequences for the EUV flux distri- 
bution from such stars (e.g. Sellmaier et al. 1996; Schaerer 
& de Koter 1997). 

We have therefore used the grid of theoretical stellar 
flux distributions from Schaerer & de Koter (1997). In con- 
trast with the Kurucz (1991) O star flux distributions, those 
of Schaerer & de Koter (1997) (i) differ greatly below the 
He II A228 edge, and (ii) show little metallicity dependence in 
the important AA228-912 region. Indeed, Kurucz (1991) and 
Schaerer & de Koter (1997) flux distributions show only mi- 
nor differences at low metallicities in this region. Stasihska & 
Schaerer (1997) have demonstrated the improved agreement 
between theoretical predictions and observations of H II re- 
gions for Schaerer & de Koter (1997) fluxes compared to 
those from Kurucz (1991). 

The Schaerer & de Koter (1997) have an additional ad- 
vantage over Kurucz in that they are coupled to the Geneva 
stellar evolution code, so that each flux distribution relates 
to a particular age, luminosity and mass. We shall there- 
fore use the main-sequence (A2-F2) models from Schaerer 
& de Koter (1997), parameters of which are listed in Ta- 
ble ^. These were calculated with a metallicity of 1/5 Zq, 
appropriate to the metal content of NGC 5253. Two tem- 
perature entries provide equivalent O spectral types from 
these models: Vacca et al. (1996) constructed an empirical 
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T ett — calibration based on observational and theoretical re- 
sults available at the time. An updated OB calibration was 
provided by Crowther (1997) including more recent results, 
and omitting the linear T e ff-spectral type relationship that 
was assumed by Vacca et al. (1996). 



4.4 Extinction Effects 

We cannot assign a single extinction to the central region 
of NGC 5253, because the extinction towards the center is 
clearly very patchy and irregular (as can be seen in the dust 
lane and in the Ha/Hf3 ratio map of Calzetti et al 1997). 
Observers have found extinction values which range from 
A v < 1 to A v = 35 mag depending on the wavelength, aper- 
ture, and exact location of the observation. The center of 
the galaxy is an association of many small regions of greatly 
differing extinctions; we have here to determine reasonable 
values of extinction for the hot nucleus and the cooler core. 

Calzetti et al. (1997) find that the true nucleus, which 
we have identified with the [S iv] source, is quite opaque in 
the visible; comparing the Ha flux to the 2-cm radio they 
found that the nuclear star cluster is embedded in a cloud 
of total ^4„=9-35 magnitudes. The Brackett a and 7 lines of 
Kawara et al (1989) give A v of 7-12 magnitudes in an area 
of 10x20 arcsec. The extinction to the rest of the core region 
is much lower, such that the total Ha compared to the 6-cm 
flux in the inner 38" gives A v of about 1. This agrees with 
the model that the youngest stars in NGC 5253 are in the 
nucleus; it is expected that the youngest stars will be the 
more highly obscured as they have not had time to disperse 
the dense molecular cloud in which they formed. 

For this work we shall use the Kawara et al. (1989) 
Br7/Bra ratio since the region sampled by Kawara et al. 
is comparable to the ISO aperture. The relation between 
the extinction in the infrared and A v is also not well es- 
tablished. Nevertheless, one magnitude of 9.7^im extinction 
is estimated to be between 15 and 30 magnitudes A v . Al- 
though the extinction in the infrared is much lower than 
in the optical it can still have significant effects on the line 
ratios. Extinction curves in the infrared differ from source 
to source, but in broad terms the [S iv] and [Ar ill] lines 
will be most affected, and similarly affected, by extinction. 
From the Kawara et al. Br7/Bra ratio we find A v =7.7 mag, 
^4k=0.7 mag and j4m=0.2 mag using the Howarth (1983) 
extinction curve. Consequently, the extinction at 9.7/xm, is 
0.25-0.5 mag, and from the curve of Draine (1989) there is 
0.2-0.4 mag extinction at [Arm] and [Siv], 0.1-0.2 mag at 
[S in], 0.08-0.17 mag at [Nell], and 0.06-0.13 mag at [Nem]. 



5 RESULTS OF PHOTO-IONIZATION 
MODELLING 

Following the approach described above, we now proceed to 
determine the massive stellar content in both the starburst 
nucleus and core of NGC 5253. Fig ^ shows how selected IR 
line ratios are affected by changes in temperature and ion- 
ization parameter. For a given density, the dependance of 
line ratios on Q translates into a sensitive function of the 
input ionization, and thus on the stellar temperature. It is 
apparent that [S iv] emission is predicted at high Q and T e s , 



Table 3. (a) Results from photo-ionization models of individ- 
ual O stars (Schaerer &: de Koter 1997) that are consistent with 
selected diagnostic IR line ratios for the nucleus of NGC 5253. 
Preferred solutions are indicated by • 





log Q 


Ar III 


Siv 


No in 


Siv 


Nem 


Siv 




kK 




Siv 


Bra 


Siv 


Ne 11 


Ne 11 


Sin 




Observed 


0.09 


~25 


<1.4 


>39 


>46? 


>1.4 




36 


9.0 


0.08 


11 


0.8 


37 


29 


1.8 




38 


8.25 


0.10 


11 


1.1 


26 


30 


1.8 




38 


8.5 


0.07 


15 


0.9 


55 


48 


2.6 


• 


38 


8.75 


0.05 


20 


0.7 


115 


76 


3.5 




42 


7.75 


0.15 


7 


1.4 


13 


19 


1.0 




42 


8.25 


0.09 


15 


0.7 


56 


48 


2.5 


• 


42 


8.75 


0.04 


21 


0.5 


155 


83 


4.2 




46 


7.75 


0.13 


8 


1.3 


18 


22 


1.2 




46 


8.25 


0.06 


15 


0.7 


71 


50 


2.6 


• 


46 


8.75 


0.04 


23 


0.5 


182 


98 


4.7 




50 


7.75 


0.13 


9 


1.2 


20 


24 


1.3 




50 


8.25 


0.06 


16 


0.7 


78 


54 


2.8 


• 


50 


8.75 


0.04 


23 


0.5 


190 


100 


4.8 





Table 3. (b) Identical to Table 3a, except for the core of 
NGC 5253. 



Tcff log Q 


Arm 


Ne H 


Si II 


Sni 


Siv 


kK 


Neil 


Bra 


Nell 


Neil 


Neil 


Observed 


0.08-0.3 


~3.3 


1.2-2.2 


<2 


<0.35 


34 7.25 


0.08 


3.2 


0.8 


1.3 


0.01 


34 8.0 


0.13 


2.6 


0.5 


2.3 


0.09 


34 8.75 


0.29 


1.8 


0.1 


3.9 


0.62 


35 7.0 


0.13 


3.0 


1.0 


1.3 


0.03 


35 7.5 


0.17 


2.5 


0.9 


2.1 


0.11 


35 8.0 


0.31 


1.7 


0.8 


3.5 


0.47 


36 6.0 


0.11 


3.6 


1.6 


0.4 


0.00 


36 6.5 


0.15 


3.3 


1.3 


0.8 


0.02 


36 7.0 


0.21 


2.7 


1.2 


1.5 


0.09 


37 6.25 


0.17 


3.5 


1.5 


0.7 


0.02 



while [Nell] originates from lower Q and T" e ff regions. (Al- 
most identical contours are obtained using the solar metal- 
licity theoretical flux distributions of Schaerer & de Koter 
(1997)). We shall now exploit these dependances to derive 
characteristic stellar temperatures and nebular densities for 
the two regions under consideration. 



5.1 [Siv] emitting central nucleus 

In Section ^ we argued that the [S iv] emission is concen- 
trated in the nucleus region of NGC 5253 observed by Beck 
et al. (1996). Our principal nucleus diagnostics are therefore: 
(i) [Arm] /[Siv] from Beck et al. (1996); (ii) [Siv]/[Neii] 
from Beck et al. (1996); (iii) [Siv]/Bra, given that ~30% 
of the total ISO Bra flux is emitted in the nucleus (Sec- 
tion 3.5); (iv) an upper limit on [Ne in] /[S iv] from ISO, 
since [S iv] is formed exclusively in the core; (v) a lower 
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Figure 4. Sensitivity of selected IR line ratios to ionization parameter (Q) and stellar temperature (T e ff). Circles indicate observed line 
ratios in the nucleus (filled-in) and core (open), while squares indicate ratios for the combined core-nucleus region. Some observed ratios 
lie within a particular range of values (connected by solid lines) while others are merely limits (indicated by arrows). 



limit on [Siv]/[Siii] for similar reasons. [Nem]/[Neii] rep- 
resents a secondary diagnostic since it relies on assumptions 
regarding the formation of [Nelll] in the nucleus. 

In Fig. filled-in circles indicate the position of these 
IR line ratios in the nucleus. Each diagnostic supports 
a high ionization parameter for the nucleus. For exam- 
ple, [Ar m]/[S iv] requires solutions in the range (T e g/kK, 
log Q)=(36, 9) to (50, 8) from Fig. |(f). The observed 
[Siv]/Bra ratio suggests even higher ionization parameters 
from Fig. ^(a). 

Unfortunately, despite our best efforts, no single solu- 
tion exists for which all observed line ratios are satisfied. As 
previous investigators have also found, it is extremely diffi- 
cult to break the degeneracy between ionization parameter 
and hardness of the ionizing spectrum from modelling alone. 
At most two or three are in agreement, suggesting errors in 
our assumptions relating to line formation regions, photo- 
ionization modelling or flux distributions. Table |^ compares 
selected IR diagnostic ratios at particular temperatures and 



ionization parameters with observed ratios for the nucleus. 
We find that optimum agreement is obtained for a narrow 
'corridor' ranging in (T e ff/kK, log Q) from (38, 8.5) to (50, 
8.25). Each of these solutions reproduce the observed neb- 
ular properties equally well. Consequently, the nucleus is 
dense and almost certaintly contains very early O-type stars. 
The minimum nucleus stellar temperature corresponds to 
an 07 V star and higher temperatures to earlier types. As 
stated above, the ionization requires 1000 equivalent 07 V 
stars; if the actual stars are younger, fewer will be needed. 

Our principal diagnostic ratios have the basic problem 
that they combine ions of two different elements and may 
therefore be in error if the elemental abundances determined 
by Kobulnicky et al. (1997) are incorrect. This is why ratios 
of lines from different ions of the same element are usually 
to be preferred. However, in NGC 5253 the two emission 
regions appear to be in such extreme states of temperature 
and ionization parameter that the contribution of one ionic 
line completely dominates. 
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5.2 [Nell] emitting core 

We now discuss constraints on the stellar content and den- 
sity of the [Ne n] emitting core regi on o f NGC 5253. Since the 
[Aril] line flux is uncertain (Sect. 2.2) we are unable to use 
the [Arii]/[Nen] ratio (<0.35) to constrain the core proper- 
ties. All solutions lie in the range 0.15<[Ar n]/[Ne n]<0.36, 
assuming a Ne abundance from Kobulnicky et al. (1997) and 
Ar abundance of 1/4 solar, 

Consequently, our choice of diagnostics for the core re- 
gion are: (i) [Nell]/Bra from ISO, adjusted for the presence 
of Bra in the nucleus; (ii) [Arm] /[Nell], adjusted for the 
presence of [Arm] in the nucleus and the maximum range 
in line fluxes; (iii) [Sin] 34.8/im/[Ne n], using ISO measure- 
ments from Genzel et al. (1998), assuming the Si abundance 
from Kobulnicky et al. (1997). Note that these fluxes were 
obtained from different sized apertures, so this ratio should 
be treated with caution; (iv) An upper limit on [S ill] /[Nell]. 
(In reality, [S ill] is anticipated to form in both the nucleus 
and core regions, analogous to [Arm]). We must wait for 
mid-IR instruments such as Michelle for observational ver- 
ification of these predictions. The [S iv] flux originating in 
the core is formally < lxlO -16 Wm -2 , obtained by sub- 
tracting the nucleus [S iv] flux of Beck et al. (1996) from the 
maximum possible ISO flux. However considering sources of 
uncertainty, a strict upper limit of < 4xl0 16 Wm -2 is cho- 
sen, so that we may additionally use an upper limit on the 
[Siv]/[Nen] ratio. 

Observed IR line ratios for the core are indicated by 
open circles in Fig. ^ while Table ^|(b) summarises selected 
(T e s,Q) combinations. Once again, no single solution exists 
for which all our diagnostics are reproduced. Nevertheless, 
a very low stellar temperature of 34-36kK is favoured for 
the ionizing O stars, with the ionization parameter poorly 
constrained. 

In order to produce sufficient Lyman continuum ioniz- 
ing photons for the observed Bra flux in the core requires 
2 500 equivalent 07 V stars (Sect. 3.2) - significantly more 
if the stars are much cooler, as appears to be the case: 
T e ff ~35kK corresponds to a stellar type of 08 or 09 (Ta- 
ble g) . Table |(b) appears to permit ionization parameters as 
low as log Q~6. However, for this Q-value, the entire volume 
of the core would only just be sufficient to contain the H II 
regions associated with 2,500 07 V equivalents. Therefore 
a higher ionization parameter is favoured, albeit probably 
lower than the nucleus of NGC 5253. Stasihska & Leitherer 
(1996) obtained a fairly tight spread of ionization parameter 
for young starburst regions. 

Combining the above material gives a picture of a hot, 
dense nucleus with T e ff > 38kK and log Q ~ 8.25, contain- 
ing 1 000 07 V "equivalents" , plus a cooler core of T e ff~35kK 
and log Q <8, with 2 500 07V "equivalents". It should be 
noted that the core is cool only by comparison with the nu- 
cleus; in a normal star forming system the low ionization 
lines of [Nell] and [Aril] dominate the mid-IR spectrum, 
which is not the case in NGC 5253. 

Our results thus far have relied on the validity of 
the photo-ionization model CLOUDY. Although this code is 
widely employed, and has been verified against a suite of 
other photo-ionization codes (Ferland et al. 1995), we have 
performed additional calculations using the shock photo- 
ionization code mappings (Sutherland & Dopita 1993). The 



mappings results agree with the overall picture of the hot 
nucleus in the cooler core, but find rather higher ionization 
parameters and stellar temperatures for both regions. We 
believe this discrepancy can be principally attributed to the 
use of different O star flux distributions (Kurucz in the case 
of mappings) and may serve as a reminder - and measure 
- of the disagreements currently inherent in the use of dif- 
ferent photo-ionization models. To verify our results further 
we have performed additional calculations, using theoretical 
Wolf-Rayet flux distributions. 



6 CONTRIBUTION OF WOLF-RAYET STARS 

Up to now we have assumed that O stars are solely respon- 
sible for the Lyman ionizing flux from the entire galactic 
core. However, NGC 5253 is a Wolf-Rayet galaxy, with both 
WN-type and WC-type stars present (Schaerer et al. 1997). 
Schaerer et al. estimated a population of <25 WN stars and 
10 WC stars in the nucleus sampled by the ground-based 
IR data sets of Beck et al. (1996), with an additional 40 
WN stars and 13 WC stars located 3" (60pc) to the south 
(also within the ISO beam). Schmutz, Leitherer & Gruen- 
wald (1992) have indicated that Wolf-Rayet stars produce a 
harder flux distribution than O-type stars, with significant 
energy emitted shortward of A=228A, except for those stars 
with the strongest winds. Is it possible that the WR stars of 
NGC 5253, though few in number, contribute significantly to 
the high excitation nebular lines (e.g. [Siv]), and invalidate 
our earlier conclusions? 



6.1 Wolf-Rayet flux distributions 
line blanketing? 



what role for 



The widely used WR energy distributions of Schmutz et al. 
(1992) employ non line-blanketed, pure helium, model atmo- 
spheres. Models including light metals (e.g. CNO elements) 
are also in use (e.g. Crowther, Smith & Hillier 1995). Line 
blanketing by heavy elements (especially Fe) is known to 
dramatically affect the emergent extreme UV flux distribu- 
tions (see e.g. Crowther, Bohannan & Pasquali 1998), but 
only test model calculations with heavy metals have been 
published to date (Schmutz 1997; Hillier & Miller 1998), 
and large grids not yet available. We have therefore cal- 
culated representative line-blanketed models for late WN- 
type ( WNL) , strong lined early WN-type (WNE-s) and early 
WC-type (WCE) stars at low metallicities (~0.2 Zq) using 
the Hillier & Miller (1998) approach. We have omitted cal- 
culations for late WC stars since they are not considered to 
exist in low metallicity environments. The reader is referred 
to Hillier & Miller (1998) for details of the technique used, 
except that our calculations include hydrogen, helium, car- 
bon, nitrogen (for WN stars), oxygen (for WC stars) and 
iron. 

Fig. g presents a comparison between the EUV flux dis- 
tributions of these WR models with the Schaerer & de Koter 
(1997) models of low metallicity dwarf O stars. WNE-s mod- 
els incorporating line blanketing show flux distributions in 
the 228-912A range which are fairly similar to the 50kK 
O star model of Schaerer & de Koter (1997) at low metal- 
licity. The strong, line-blanketed WR stellar wind hinders 
emission below A=228A, in contrast with the most recent 
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Figure 5. Comparison between latest Lyman continuum flux dis- 
tributions from O stars (Schaerer & de Koter 1997) at low metal- 
licity (0.2Zq) with recent line blanketed model atmospheres for 
Wolf-Rayet stars (following Hillier & Miller 1998). It is apparent 
that strong-lined WCE and WNE stars show similar EUV flux 
distributions to hot O stars, except that negligible flux is emitted 
at high energies (3—4 Ryd), while solely the WNE-w model shows 
a strong flux above 4 Ryd (A < 228A). Ionization edges relevant 
to this study are indicated. 

recent star models. WCE models are similar except that 
negligible flux is emitted A <300A, due to the blanketing by 
carbon and oxygen. WNL stars show a much softer energy 
distribution. 

Our results appear in conflict with the much harder 
fluxes predicted by certain unblanketed models of Schmutz 
et al. (1992), particularly below the He + edge. However, 
Schmutz et al. correctly stress the importance of stellar 



Table 4. Comparison between IR line ratios predicted in photo- 
ionization modelling of NGC 5253 using line blanketed O-type 
models (Schaerer & de Koter 1997) and WR distributions (fol- 
lowing Hillier & Miller 1998) at O.2Z . The O star T cff -spectral 
type calibration of Crowther (1997) is adopted. 
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wind density, such that emission at A <228A relies on the 
WR wind being relatively transparent. Denser winds de- 
stroy photons beyond this edge. Does this remain valid 
for metal-line blanketed models? In Fig. ^ we include the 
predicted EUV spectrum for a line blanketed model of 
the extremely weak-lined, early-type WN (WNE-w) star 
HD 104994 (WR46, WN3) using parameters from Crowther 
et al. (1995), except for a low Fe-abundance of 0.2Zq. In 
agreement with Schmutz et al. the flux distribution for this 
star is in dramatic contrast with other WR (and O star) line 
blanketed models in that it displays a very hard spectrum 
above the He + edge. This is be expected, since a subset of 
WR stars do emit a significant flux shortward at A <228A, 
as evidenced from nebular He II A4686 emission (e.g. Dopita 
et al. 1990; Garnett et al. 1991). 

To test for the influence of WR stars we have carried 
out CLOUDY calculations similar to those discussed above. In 
Table ^ predictions for selected IR line ratios at log Q—6.5, 
7.5 and 8.5 using our WR energy distributions are contrasted 
with those from a variety of O star models. These results are 
now discussed below. 

For all realistic ionization parameters, WNL (specifi- 
cally WN8) stars have negligible contribution to the high ex- 
citation nebular lines observed in the nucleus of NGC 5253, 
as may be expected from Fig. ^| WNL stars do, however, 
make contributions to the low excitation features ([Nell], 
[Sin]) in the core of NGC 5253 that are comparable to late 
O dwarfs. We can safely conclude that the presence of 40-65 
WNL stars in NGC 5253 (according to Schaerer et al. 1997) 
will not significantly affect our earlier results based solely on 
O stars. 
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Figure 6. Predictions of flux distributions for a burst of star 
formation at an age of 2— lOMyr from Schaerer & Vacca (1998) 
at Z=0.2Zq, using the Meynet et al. (1994) evolutionary tracks 
coupled with O star flux distributions from Schaerer &; de Kotcr 
(1997), plus Schmutz et al. (1992) WR models, present between 
2.8— 5.2Myr (causing the hard EUV flux distribution). Ionization 
edges are indicated. 



The similarity in the EUV flux distributions of WNE- 
s and WCE models to early O dwarfs is reflected in their 
comparable diagnostic line ratios of neon and sulphur. Con- 
sequently, the presence of ~23 WC stars in NGC 5253 will 
not influence our results for either region. 

In general, the presence of WR stars in young starburst 
regions do not affect results from exclusively O stars (how- 
ever see below). WNL flux distributions do not differ sig- 
nificantly from very late O stars, while WNE-s and WCE 
energy distributions closely mimic early-type O stars. As we 
shall see in Sect. [t], this is in sharp contrast with recent evo- 
lutionary synthesis models (Schaerer & Vacca 1998) which 
utilise the WR models of Schmutz et al. (1992) showing 
strong emission at A <228A. 



6.2 Wolf-Rayet stars as the source of [O iv] 
emission in starbursts? 

Lutz et al. (1998) have recently proposed the identification 
of [0 iv] 25.91/^m emission in several young starburst galax- 
ies, including NGC 5253. They suggest that WR stars may 
provide the hard ionizing flux necessary for such a highly 
excited feature. Our ISO spectrum suggests an upper limit 
for the [O iv] line of 0.1 x KT 16 W m~ 2 or [O lv]/Bra<0.07, 
reasonably assuming any emission is restricted to the nu- 
clear region. (Lutz et al. report 0.65 xlO -16 W m -2 based 
on higher S/N data sets, implying [O lv]/Bra~0.35). Our 
observations indicate that [Nev] 14.3/im is not present in 
NGC 5253. 

Line blanketing softens the flux from the majority of 
early-type WR stars, such that negligible nebular [O iv] 
25.9/^m emission is predicted from photoionization models 
of such stars, in contrast with early O stars (Table |^). This 
is the precise opposite of predictions using earlier O and WR 
models (Kurucz 1991; Schmutz et al. 1992)! Within the high 
ionization parameter nucleus of NGC 5253, the low metal- 
licity grids of Schaerer & de Koter (1997) predict that ob- 
servable nebular O iv could be produced by early O dwarfs: 
[Oiv]/Bra~0.5 for 03 V stars at log Q ~8.5. 

However, certain WR stars emit large numbers of pho- 
tons with energies above the Hell A228 edge, necessary for 
[O iv] emission. We include results of photo-ionization mod- 
elling of the WNE-w flux distribution in Table From 
the table, [O iv] 25.9/^m emission is predicted at all ioniza- 
tion parameters, two orders of magnitude greater than 03 V 
stars at parameters appropriate for the nucleus of NGC 5253! 
Other high ionization IR nebular lines are predicted, includ- 
ing [Nev] 14.3^m ([Nev]/Bra~6 at log Q=8.5). 

Since such emission features are not prominent in our 
ISO observations, we can exclude the possibility of a sig- 
nificant number of such stars being present in NGC 5253. 
This is expected, since WNE-w stars are very rare in our 
Galaxy and the Magellanic Clouds. They have unusually 
low emission line fluxesjj]. Consequently, their identification 
would prove difficult in external galaxies from optical spec- 
troscopic techniques. 

Our results suggest that WR stars will only play a major 
role in the overall energy distribution of a young starburst 
region if they are very hot, and possess stellar winds that are 
relatively transparent. In the nearby universe such stars are 
rare, so it is unlikely that they represent the 'Warmers' re- 
sponsible for the high IR excitation features observed in Ac- 
tive Galactic Nuclei (Terlevich & Melnick 1985). However, 
metallicity also complicates matters yet further, affecting 
both the degree of EUV-line blanketing and the strength 
of the stellar wind (although a metallicity dependence of 
mass-loss in WR stars has yet to be established observa- 
tionally). If such a dependence is identified, weaker stellar 
winds would yield harder EUV flux distributions. Such a 
result would have an important effect on the possibility of 
WR stars providing the hard ionizing flux in very metal- 
poor starbursts, such as I Zw 18 (de Mello et al. 1998). 



t WR46 has a Hell A4686 emission line flux of 1.2xl0 35 ergs -1 
which is a factor of 4 lower than typical WNE stars (Schaerer & 
Vacca 1998) 
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Table 5. (a) Results from photo-ionization models of evolution- 
ary synthesis calculations for an instantaneous burst at 0.2Zq 
(Schaerer & Vacca 1998) that are consistent with selected diag- 
nostic IR line ratios for the nucleus of NGC 5253. Our preferred 
solutions (indicated by •) are prior to the hot WR phase (2.9— 
4.7 Myr) because of the uncertain WR flux distributions assumed 
in the models. 



Age 


log Q 


Ar III 


Siv 


Ncin 


Siv 


Ncin 


Siv 


Myr 




Siv 


Bra 


Siv 


Neil 


Neil 


Sm 


Observed 


0.09 


~25 


<1.4 


>39 


>46? 


>1.4 


2.0 


8.0 


0.10 


10 


1.0 


11 


30 


1.7 


2.0 


8.25 


0.07 


1 1 


0.7 


60 


45 


2.4 • 


2.0 


8.5 


0.08 


18 


0.6 


102 


63 


3.3 


2.5 


8.0 


0.11 


9 


1.1 


21 


24 


1.3 


2.5 


8.25 


0.08 


12 


0.8 


43 


35 


2.0 • 


2.5 


8.5 


0.06 


16 


0.7 


79 


52 


2.8 


3.0 


8.0 


0.09 


10 


1.0 


34 


32 


1.7 


3.0 


8.25 


0.07 


14 


0.7 


66 


48 


2.5 • 


3.0 


8.5 


0.05 


18 


0.6 


117 


71 


3.5 


4.0 


8.0 


0.10 


10 


1.0 


30 


30 


1.6 


4.0 


8.25 


0.07 


13 


0.8 


59 


45 


2.3 


4.0 


8.5 


0.05 


16 


0.6 


100 


63 


3.0 


4.5 


8.25 


0.09 


10 


0.9 


35 


32 


1.7 


4.5 


8.5 


0.08 


11 


0.8 


68 


33 


2.0 


4.5 


8.75 


0.05 


17 


0.6 


112 


68 


3.4 



Table 5. (b) Identical to Table 5a, except for the core of 
NGC 5253. 



Age 


logQ 


Arm 


Neil 


Si H 


Sm 


Siv 


Myr 




Neil 


Bra 


Neil 


Neil 


Neil 


Observed 


0.08-0.3 


~3.3 


1.2-2.2 


<2 


<0.35 


4.6 


6.25 


0.18 


3.5 


1.6 


1.0 


0.02 


4.7 


6.0 


0.12 


4.2 


1.7 


0.4 


0.01 


4.7 


6.5 


0.18 


3.2 


1.5 


0.9 


0.03 


4.7 


7.0 


0.29 


2.4 


1.4 


1.9 


0.19 


5.0 


6.75 


0.10 


3.3 


1.1 


0.9 


0.02 


5.0 


7.5 


0.14 


2.9 


0.7 


1.8 


0.13 


5.0 


8.25 


0.27 


2.1 


0.4 


3.1 


0.83 


5.2 


7.5 


0.08 


3.6 


0.5 


1.3 


0.05 


5.2 


8.0 


0.11 


3.7 


0.3 


1.7 


0.15 


5.2 


8.5 


0.19 


2.7 


0.2 


2.4 


0.55 



More definitive answers await the calculation of a grid of line 
blanketed WR models for various combinations of mass-loss 
rates and metallicity. Crowther (1999) discusses this aspect 
in greater detail. 



7 PHOTO-IONIZATION PREDICTIONS FOR 
A BURST OF STAR FORMATION 

Our photo-ionization modelling has used a single O or WR 
star flux distribution, adjusted to differing ionization param- 
eters. In reality, the nucleus and core of NGC 5253 will con- 



tain stars spanning an enormous range of masses and tem- 
peratures, even if they were formed in a single burst of star 
formation. So, how valid is our assumption that the most 
massive (and consequently the highest temperature) stars 
dominate the IR nebular line fluxes? To investigate this, we 
have kindly been provided with evolutionary synthesis mod- 
els (v 2.32) from Schaerer & Vacca (1998). These utilise the 
low metallicity (0.2Zq) Schaerer & de Koter (1997) flux dis- 
tributions for O-type stars, a Salpeter IMF, a mass range of 
0.8-120Mq and were calculated for differing ages accord- 
ing to the Meynet et al. (1994) evolutionary tracks after an 
instant burst of star formation. 

Fig. ^ presents EUV flux distributions produced by 
these models for different ages. Excess high energy photons 
at A <228A are produced by WN and WC stars between 
2.8-5.2 Myr, since Schaerer & Vacca (1998) utilised the low 
wind density Schmutz et al. (1992) WR model distributions. 
From Section |^ such EUV excesses may be error. Neverthe- 
less, detailed comparisons may be made with our predictions 
based on single O star models, treating specific results dur- 
ing the WR phase with caution. We have generated a grid of 
photo-ionization models for ages in the range 2-10 Myr and 
calculated fluxes for a fixed density of 10 cm -3 , with the 
total luminosity running from 10 1 , 10 2 , 10 n L©, which 
covers a range in ionization parameter similar to our ear- 
lier calculations. The sensitivity of IR diagnostic line ratios 
to ionization parameter and age (in Myr) are presented in 
Fig. 0, where we again include observed line ratios in the nu- 
cleus (filled circles), core (open circles) and combined core- 
nucleus region (open squares). 

In Table ^(a) we present selected age and ionization 
parameter solutions that are in reasonable agreement with 
our IR diagnostics in the nucleus. The degree of consistency 
between alternative diagnostics is comparable to those from 
individual O star flux distributions (Table §[a]). Our ear- 
lier results for the nuclear region ionization parameter are 
supported, with log Q~8-8.5. The precise age is difficult 
to tightly constrain, but ages in the range 2-4 Myr are in 
reasonable agreement with the observed line ratios. Recall- 
ing that the synthesis calculations between 2.9-4.7 Myr are 
strongly affected by the presence of hot WR stars, our pre- 
ferred solution is 2-3 Myr. 

For the core, Table |B|(b) presents results of our photo- 
ionization models that are in reasonable agreement with our 
chosen mid-IR diagnostics, though the level of agreement is 
no better than that obtained from individual O star flux 
distributions (Table ^[b]). Once again, we find that the core 
ionization parameter remains uncertain, with log Q=6.5 to 
8 providing reasonable agreement for a fairly narrow range 
of ages, between 4.7 to 5.2 Myr. The only Wolf-Rayet stars 
predicted during this age range are late WN stars, with a 
fairly soft EUV distribution, so that the evolutionary syn- 
thesis models should be appropriate. A significantly older 
population (5-10 Myr) would require a very low ionization 
parameter (log Q~6). As previously indicated, the entire 
volume of the core would only just be sufficient to contain 
the Hn regions associated with 2,500 07 V equivalents for 
this Q value, so we can exclude such a population. 

Overall, our results from earlier calculations are broadly 
supported here, based on more sophisticated flux distribu- 
tions. Consequently, an important result from this work is 
that meaningful results may be obtained for very young star- 
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burst regions based on individual O star flux distributions 
alone. 



8 COMPARISON WITH PREVIOUS 
DETERMINATIONS 

We now compare our present results for the hot, massive 
stellar content in NGC 5253 with those obtained previously 
from ground (Aitken et al. 1982; Beck et al. 1996) and space- 
based (Lutz et al. 1996) IR observations, as summarised in 
Table @. (We are unaware of any studies based on optical 
data sets). 

All previous mid-IR analyses relied on: (i) O-star flux 
distributions which did not account for non-LTE effects or 
stellar winds, or in many cases the low metallicity of the 
galaxy; (ii) inappropriate line ratios (given the results from 
ground- and space- based observations), and (iii) assumed 
elemental abundances (often solar composition) and ioniza- 
tion parameter. Consequently, previous approaches may lead 



to imprecise results, although all studies support the pres- 
ence of hot, young O-type stars in its central starburst, with 
characteristic temperatures ranging from 40kK to >50kK. 

Our comparison with previous results well illustrates 
the effect of different O-star flux distributions and elemen- 
tal abundances, even with an identical photo-ionization code 
and IR line diagnostics being employed. For example, Lutz 
et al. (1996) obtained a characteristic O star temperature 
of 48.5kK from the observed [Neni]/[Nell] ratio since they 
adopted solar Kurucz (1991) models, solar elemental abun- 
dances and log Q—S. From our analysis, neglecting the spa- 
tial information for line formation in NGC 5253 and adopt- 
ing log Q—S implies a stellar temperature of 36kK for 
the combined core-nucleus region using [Ne m]/[Ne n] from 
Fig. ^(e). Alternatively, a stellar temperature of 39kK would 
result from the observed [Siv]/[Siii] ratio with an identical 
ionization parameter (Fig|^[b]). 

Having established constraints on the stellar content of 
the nucleus and core of NGC 5253 we are able to compare 
ages obtained from the derived stellar temperatures with 
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Table 6. Summary of comparison between stellar temperatures and ionization parameters derived for NGC5253 using IR diagnostics. 
Theoretical O star fluxes are from Hummer & Mihalas (1970, HM70), Kurucz (1991, K91) and Schaerer & de Koter (1997, SdK97), while 
photo-ionization models employed are MAPPINGS II (Sutherland & Dopita 1993, SD93) and CLOUDY (Ferland 1996, F96). 



Study 


IR Diagnostics 


- O star Fluxes - 
Ref. Z Q 


Photo- 
Code 


Abund. 

Z & 


kK 


logQ 


Age 
Myr 


Notes 


Aitken et al. 1982 


[Siv] 


HM70 


1.0 




1.0 


>50 








Beck et al. 1996 


[Siv], [Arm], [Neil] 


K91 


0.1-1.0 


SD93 


1.0 


40-45 








Lutz et al. 1996 


[Neii-m] 


K91 


1.0 


F96 


1.0 


48.5 


8.0 




Q assumed 


This work (nucleus) 


[Siv], [Arm], [No in], Bra 


SdK97 


0.2 


F96 


0.25 


38-50 


~8.25 


2-3 




This work (core) 


[Neil], [Arm], [S III], Bra 


SdK97 


0.2 


F96 


0.25 


~35 


<8 


4.7-5.2 





those from evolutionary synthesis calculations. A minimum 
characteristic stellar temperature of >38kK for the nucleus 
implies a minimum characteristic mass of ~35M© . Using the 
result of Kunze et al. (1996) that the upper mass cut-off of a 
cluster is about 5Mq greater than the single-star equivalent 
mass gives a lower limit to the mass cutoff of 40Mq for 
the hot nucleus, with an age no greater than around 3 Myr 
(Meynet et al. 1994). This is in reasonable agreement with 
the evolutionary synthesis calculations for log Q~8, and is 
typical of recent age estimates for the central nucleus of 
NGC5253 from the WR population (Schaerer et al. 1997). 
(For T c ff ~46kK the implied hottest stellar type is 04 or 05 
with mass ~60A/© at a corresponding age of ~2.5Myr). 

For the core region, a stellar temperature of ~35kK 
corresponds to the hottest spectral type being equal to or 
later than 08V or 09V (Table ^), and a equivalent mass 
of ~25M0. This gives an upper mass cutoff of 30M@, indi- 
cating an age of around 4-5 Myr, in good agreement with 
that implied from our synthesis calculations. However, we 
note that our determination is based on the entire core 
of NGC 5253, spanning several hundred parsec. We do not 
claim that the massive stellar population in the core of 
NGC 5253 was formed coevally ~5 Myr ago. Our results are 
almost certainly biased towards the youngest, most massive 
clusters exterior to the central nucleus. In reality, the con- 
tent of the core will span a large range of ages. The weakness 
of [Ar n] emission does appear to rule out a particularly old 
population. Vanzi & Rieke (1997) concluded from near-IR 
spectroscopy, that NGC 5253 is young (~8Myr) and con- 
tains high temperature stars (>40kK) since He I 1.70/im is 
present, [Fell] emission (from SN remnants) is very weak, 
as are the CO absorption bands (produced by red giants). 
Davies et al. (1998) estimated a greater age of ~10-12Myr 
for regions ~7arcsec from the central nucleus. 

Again, neglecting the spatial origin of low and high 
excitation lines and assuming a uniform ionization param- 
eter of log Q=8, an age of ~4.8 Myr is obtained from 
[Ne ill] /[Neil] following Fig. 0(e). Alternatively, ages of ei- 
ther 2.9 or 4.4 Myr are implied from the observed [S iv]/[S ill] 
ratio (Fig. [}b]). 

To summarize, many models and techniques may be ap- 
plied to IR observations of this galaxy. While they differ in 
details they agree with each another in the general picture 
of NGC 5253. This galaxy has a very young starburst in the 
nucleus and a slightly older one within the core. It is remark- 
able for the youth of the starburst, probably no more than 
3 Myr, and for the very massive stars it formed. (For stars to 



reach the WC phase in such a low metallicity environment 
requires very massive progenitors of 6O-12OM0 (Maeder & 
Meynet 1994).) 



9 DISCUSSION AND CONCLUSIONS 

NGC 5253 is widely acknowledged to be an excellent, in some 
respects unequalled, laboratory for the study of intense star- 
bursts. This is because it is relatively close, contains a very 
young starburst, and because its central region is so varied, 
containing many sub-regions and clusters of different ages 
and stellar content. (The patchiness of the central region of 
NGC 5253 is probably related to its extreme youth; there has 
not been time for stars to migrate away from their points 
of formation and for the original structure of star forma- 
tion to be lost). NGC 5253 also has the special advantage of 
dwarf galaxies: in contrast with large spiral galaxies bursts 
are seen against a very low intra-burst star formation rate 
rather than steady star formation activity. 

It is therefore rather troubling that the properties of 
NGC 5253 at which we arrive in this paper depend so cru- 
cially on the high spatial resolution ground-based data. The 
characteristics of NGC 5253 that make it so interesting also 
make it very hard to interpret observations that could not, 
for example, distinguish between the nucleus and core. If 
we had not been able to assign the different infrared lines 
to their appropriate sources, we would have found highly 
contradictory line ratios and an imprecise effective stellar 
temperature of the ionizing O stars. Even with the exist- 
ing ground- and space- based observations, some important 
results await confirmation; e.g. What fraction of the [S ill] 
nebular flux is located in the nucleus? 

Nevertheless, it is clear from our results that the nucleus 
of NGC 5253 contains extremely young, massive, hot stars, 
with a high characteristic ionization parameter. The sur- 
rounding core contains cooler stars, with a lower ionization 
parameter, and undoubtedly contains stars spanning a wide 
range of ages. The nucleus of NGC 5253 is not unique in this 
respect; the stellar content of 30 Doradus in the LMC, for 
example, shares many of its attributes, and deserves greater 
study at mid-IR wavelengths. One important result of our 
study is that estimates of the global properties of young star- 
bursts from the effective temperatures of their most massive 
constituents appears to be a reasonable approach. 

Unfortunately it is not yet possible to readily compare 
the spatial density of O stars in the super-star nucleus of 
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NGC5253 with other giant Hn regions. Definitive answers 
await high spatial resolution IR observations, but the ion- 
ization parameter for the nucleus indicates a dense concen- 
tration. At the very least, 1000 07 V "equivalents" occupy 
a region of radius <15 parsec. For comparison, ~100 early 
O stars, with a similar total Lyman ionizing flux are lo- 
cated within a radius of 10 parsec from R136a in 30 Doradus 
(Massey & Hunter 1998; Crowther & Dessart 1998). 

Our principal result illustrates the challenge faced with 
nearly all observations of galaxies; many different sources fall 
into a single aperture. This is particularly acute for dwarf 
starburst galaxies, where there is less background of star 
formation, where the smoothing effects of shear and spiral 
arm forces are absent, and where the starburst activity is 
likely to be concentrated in one or a few clumps a few parsecs 
in size. Consequently, high spatial resolution observations of 
NGC 5253 and other nearby starbursts are urgently sought 
at near-IR and mid-IR wavelengths, and quantitative results 
for distant starbursts should, for the moment, be treated 
with caution. 
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